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Abstract - Pinning particles at random in supercooled liquids is a promising route to make 
substantial progress on the glass transition problem. Here we develop a mean-field theory by 
studying the equilibrium and non-equilibrium dynamics of the spherical p-spin model in presence 
of a fraction c of pinned spins. Our study shows the existence of two dynamic critical lines: 
one corresponding to usual Mode Coupling transitions and the other one to dynamic spinodal 
transitions. Quenches in the portion of the c — T phase diagram delimited by those two lines leads 
to aging. By extending our results to finite dimensional systems we predict non-interrupted aging 
only for quenches on the ideal glass transition line and two very different types of equilibrium 
relaxations for quenches below and above it. 



I 

C 

O 

o 



> 

00 
O 



X 

S-H 



Introduction. — The slowing down of the dynam- 
ics upon cooling is an universal property of supercooled 
liquids. In the most spectacular cases relaxation time in- 
creases by near 15 orders of magnitude in a pretty narrow 
temperature window while the structural properties of the 
ystem remain almost unchanged. Since the beginning, a 
recurrent question in the literature has been whether the 
fast growth of the relaxation time is a purely dynamic 
phenomenon or it is governed by the approaching of a 
thermodynamic singularity (the ideal glass transition) at 
a finite temperature Tk ■ The problem in resolving this is- 
sue is related to the very nature of the glass transition: the 
relaxation time increases so steeply that it is not possible 
to approach the putative transition very closely. Below a 
certain temperature, still rather far from Tk, the system 
becomes dynamically arrested — a glass — and falls out of 
equilibrium. 

A way to partially bypass these equilibration problems 
and find a method to perform stringent tests on the ideal 
glass transition and its thermodynamic character has been 
proposed in pQ. In this recent work we suggested that 
pinning a fraction c of particles of a supercooled liquid in 
an equilibrium configuration induces an ideal glass tran- 
sition, the Random Pinning Glass Transition (RPGT), 
at temperatures fairly higher than Tk- Using the Ran- 
dom First Order Transition (RFOT) theory we derived 
the phase diagram in the c — T plane, in particular a line 
of RPGTs starting from Tk, c = and ending at certain 
values Th,Ch- Actually, the procedure of pinning particles 
is a powerful probe for studying glassy systems also for 



other reasons. It has been proposed in [5J[5] and then used 
in [UEH7] as a tool for detecting the presence of medium 
range amorphous order in supercooled liquids. Moreover, 
it plays a role in studies of glassy dynamics in porous 
media. In fact, a simplified way to represent the porous 
matrix consists in freezing particles of an equilibrated liq- 
uid. The protocol leading to RPGTs correspond to the so 
called Equilibrated Mixture procedure studied in [5115]. 
In simulations, either the ones aimed at studying amor- 
phous order or the ones focusing on porous media, it was 
found that the presence of pinned particles always pro- 
duces a significative slowing down in the dynamics of the 
remaining free particles of the system [SUSIE]- This phe- 
nomenon turns out to be fully consistent with what is 
expected from our RFOT analysis pQ: the slowing down 
of the dynamics is a precursor of the RPGT. A complete 
dynamic picture including out of equilibrium dynamics is, 
however, still missing. Many related interesting questions 
remain open. For example, what happens if one quenches 
a random pinned liquid below its glass transition temper- 
ature? Is the subsequent aging dynamics the same one of 
unpinned systems? 

By analyzing the dynamics and the off-equilibrium dy- 
namics of the randomly pinned p-spin spherical model, a 
system providing a mean-field theory of glassy systems, 
wc shall answer these questions and show that indeed ag- 
ing is quite different and that this is related to the very 
special nature of the ideal glass state for RPGT. More- 
over, we shall also obtain the Mode Coupling Transition 
(MCT) phase diagram and critical properties, which dif- 
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for from the ones obtained for porous media by using pro- 
jection operators methods [TU]. A general discussion, in 
particular concerning the extension of our results beyond 
mean-field theory will be presented in the conclusion. The 
main finding of this work is that in finite dimensions ag- 
ing is always interrupted, except for quenches on the crit- 
ical RPGT line. Surprisingly, aging is interrupted even 
for quenches below the RPGT line, but very differently 
than for quenches above it. In the former case the ideal 
glass transition state is suddenly nucleated from the slowly 
evolving aging state, whereas in the latter aging smoothly 
evolves into equilibrium dynamics. These surprising re- 
sults can be tested in simulations or in experiments on 
colloids and provide new and interesting ways to ascertain 
the nature of the glass transition. 

Randomly pinned spherical p-spin model and dy- 
namical protocols. — A mean-field study of the equilib- 
rium and out of equilibrium dynamics of randomly pinned 
systems can be obtained by focusing on the spherical p- 
spin model, which is known to provide a mean-field de- 
scription of glassy dynamics. 
The Hamiltonian of the model reads 
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Fig. 1: Solutions of (|3]) for different values of the parameters c 
and T. The letter code refers the positions in the c — T phase 
diagram of FigfS] The top panel shows the development of the 
plateau in the correlation function in proximity of the dynami- 
cal transition at Td(c). The bottom panel shows that the large 
time value of the correlation function decreases without jumps 
on a path around (c h ,T h ). 
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The Si are N continuous variables satisfying the spherical 
constraint N = sf. The couplings Ji 1 ...i p are indepen- 
dent random variables extracted from a Gaussian distri- 
bution with zero mean and variance p\/2N p . 
At time t = 0, we pick a configuration C at random from 
the equilibrium measure and block a fraction c of its spins: 
Vt > 0, Si(t) = s~i for i G [l,ciV]. We focus on two differ- 
ent types of dynamics for t > 0. Equilibrium dynamics: 
the remaining free spins evolve directly from the values 
they take in the pinned configuration C: for i G [cN,N], 
Si(t = 0) = sj. In this case the system starts from (and 
remains at) equilibrium. Quench dynamics and relaxation 
toward equilibrium: the initial conditions for the remain- 
ing free spins are completely random, i.e. they assume 
the values of an infinite temperature configuration C°°: 
for i G [cN,N], Si(t = 0) = s°°. In this case we analyze 
the relaxation toward equilibrium and, when this does not 
take place, the corresponding out of equilibrium aging dy- 
namics. 

In both cases for t > the free spins evolve according to 
the Langevin equation: for i G [cN, N] 



Si(t) 



dH 

dsi 



n{t)si{t) + m(t) . 



(2) 



The Lagrange multiplier p{t) is introduced to enforce the 
spherical constraint and r)i(t) is a Gaussian white noise 
with zero mean and variance 2T. We shall be interested 
in observables averaged over the possible choices of the 
couplings, the initial configuration C and the Gaussian 
noise. 



Equilibrium dynamics. — By using the method de- 
veloped by Barrat, Burioni and Mczard [TT] to study 
the equilibrium dynamics we have obtained in the large 
N limit the equation verified by the correlation function 
C ( t > t> ) = N(i-c) E^I C Af( s »(0 s »(^))- Since the system is 
time translation invariant (it is at equilibrium), the equa- 
tion can be directly written on C(r) = C(t,t') where 
T = t - t' > 0: 



d T C{r) = -TC{t) 



f3 I dud u C(u)V'(C(T - it)) (3) 
lo 



The kernel in the integral is the derivative of V(C) 
i r„ , ri ^qy. 



(1 



2(1=3 

The solution of this equation has a shape similar to the 
one of the usual correlation function, see FigjTJ For small 
fractions of pinned spins and/or high temperatures (curves 
C2 and D in Fig[T]and corresponding spots in the phase 
diagram in Figj2]) the correlation function simply decays 
to a small value qo(c,T) oc c. Increasing c and decreasing 
the temperature, the function develops an intermediate 
plateau and shows the usual two steps relaxation of the 
correlation function in proximity of the dynamic transi- 
tion (curve Al). The plateau appears at a finite value 
(?i(c, T). Finally, along a line in the c — T plane that we 
denote either as Cd(T) or Td{c), the correlation function 
ceases to depart from qi{c,Td{c)) on any finite time, i.e. 
the system undergoes a dynamic MCT transition. This 
is the continuation in the (c, T) plane of the usual one 

that takes place at Td = J I^Efp ^r and c = 0. The 
line of dynamical transitions was already shown in [T] as 
a result of the static analysis. We plot it again in the 
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Fig. 2: Phase diagram of the p-spin model with pinned spins. 
The two dynamic transition lines meet with the thermody- 
namic transition line at (c.h,Th). Above this temperature no 
thermodynamic and dynamic transition is present in the sys- 
tem. The rectangle represents the closed path we followed to 
obtain curves of Fig[T] 

more complete phase diagran{3 of Fig 12 where we have 
also drawn the thermodynamic transition line Tr-(c) dis- 
cussed in [Tj . The large time limit of ^ gives an equation 
ong = lim-r^oo C(t): q — /3 2 (1 — q)V'(q) which for certain 
values of c and T(< Td{c)) has two solutions: qi (c, T) and 
qo{c,T). For these pairs (c, T), lim T ^oo C(r) is given by 
the highest value qi (> qo). Note that we have checked 
that replica computations lead to the same equation for 
the self-overlap q. 

The following picture results from the analysis of the dy- 
namical equation. Before the dynamic transition (points 
D, Al in Fig. 2) the system is in the paramagnetic-liquid 
phase. The small value of the correlation function at large 
time simply comes from the trivial influence of pinned 
spins towards their neighbors. At the dynamical transi- 
tion the ergodicity of the system is broken and the config- 
uration of the free spins remains trapped in a small region 
of the phase space in proximity of the pinned configura- 
tion. The large time value of the correlation function in 
this case corresponds to the self-overlap of the metastable 
state to which the pinned configuration belongs. Fig. 2 
shows that by increasing the temperature, a higher con- 
centration Cd(T) of frozen spins is required to induce the 
dynamical transition. Correspondingly, go right before 
the transition increases. Instead, the self overlap of the 
metastable state in which the dynamics remains trapped 
lowers because th e metasta ble states become broader. At 
the point T h = J ^pp-^ ^ > c h = (P ~ 2 ) 2 /P 2 at which q x 
reaches qo, the dynamical transition line ends. This cor- 
responds to a spinodal point for the metastable states, as 

lr This phase diagram is similar and in agreement with the one 
found by Ricci-Tersenghi and Scmcrjian 1 1 21 in the context of com- 
binatorial optimization problems. Note however that the transition 
taking place at T^, (c) is different from the one they discussed since 
we are focusing on non-equilibrium dynamics. 



it can be shown by the static replica computation. Above 
Th the pinning procedure is not able to induce any ergod- 
icity breaking: the high value of the large time correlation 
function is only due to trivial correlation between the free 
and the frozen spins. 

The phase diagram in Fig. 2 shows, surprisingly, that it is 
possible to follow a path in the c — T plane that connects 
the liquid and the ideal glass phase without crossing any 
phase transition, similarly to what happens for the liquid- 
gas transition. In order to illustrate this fact we discuss 
the behavior of the correlation function along the rectan- 
gular path drawn in Figf2] C(t) is singular only once, 
corresponding to the big spot in the figure. Afterwards, 
following the bottom side (A) of the rectangle the large 
time value q of the correlation function increases with c 
(see top panel of Fig. 1). On the right and the top sides 
(B and C), q decreases finally reaching, without showing 
any singularity, the low overlap solution qo of ©. At the 
same time, approaching the transition from the left (last 
part of left side and Al), the plateau forms anticipating 
the jump from qo to q\ at the dynamical transition. 
As a last issue, we address how the singular MCT behavior 
changes moving along the T<j(c) line in the phase diagram. 
The critical behavior is always the usual MCT one, called 
A2 by Gotze [13] . The only difference is in the value of the 
exponents a and b governing the approach and the depar- 
ture of the correlation function from the plateau. They 
are smaller for transitions occurring at higher concentra- 
tions of frozen spins. Accordingly, transitions with higher 
critical temperatures are characterized by more stretched 
decays towards the qo large time limit. The a and b expo- 
nent are fixed by the relation 

A _ r 2 [i -a] _ r 2 [i + b] _ 1 V"( qi (c,T d (c))) 

T[l - 2a] r[l + 2b] 2 [V"Mc,T d (c)))]3/2 ' 

For c = we recover the standard result A = 0.5, a ~ 0.395 
and b = 1. By increasing c, the values of the expo- 
nents decrease and at the point (c^T^), at which A = 1, 
they become identically zero. Correspondingly, at the end 
of the dynamical transition line, the correlation function 
shows a logarithmic decay characteristic of A 3 singular- 
ities, instead of the usual power law decay [13]. More- 
over, the relaxation time r Q diverges exponentially fol- 
lowing a generalized Vogel-Fulchcr-Tamman (VFT) law 
r a = cxp(£>|T - T h \- e ) with exponent = 1/6 Q3]. It is 
remarkable that at the only point (ch,Th) where the dy- 
namical transition does exist, since it coincides with the 
static one, mean-field theory correctly predicts a super- 
Arrhenius behavior. Finally, we point out that using the 
MCT theory of dynamical correlations [Tl] one finds that 
approaching the endpoint Th,Ch the system becomes less 
dynamically heterogeneous since the growth of X4 in time, 
given by power laws with exponents a and 6, diminishes. 
Moreover, the peak of Xi at Ch,Th diverges less rapidly. 

Non-equilibrium dynamics. — The set of equations 
describing the non-equilibrium dynamics is obtained as in 
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we have, for t > t', 



d t C(t, t')+fj,(t)C{t, t') = dt"R(t', t")V'(C(t, t"))+ (4a) 



+ dt"C(t', t")R(t, t")V"{C(t, t")) + (3C{t')V'{C(t)) 



d t R(t,t') + n(t)R(t,t') =Jdt"R(t",t')R(t,t")V"(C(t,t")) 

(4b) 

(4c) 



C(t) = p dt"R(t,t")V'(C(t")) 
Jo 



Where R(t,f) 



h i= 



is the re- 



JV(l-c) A^i=cN 6hi(t') 

sponse function (hi(t') is a magnetic field coupled to spin 
i) and the new correlation function C(t) is the correla- 
tion between the running configuration at time t and the 
pinned configuration C. By imposing the spherical con- 
straint, one obtains the evolution equation for /i(£)H 



H(t) - T = dt"R{t, t")V'(C(t, t"))+ 



(5) 



+ dt"C(t, t")R(t, t")V"(C(t, t")) + pC(t)V'(C(t)) 
Jo 

The correlation function C{t) measures how far the system 
is from reaching equilibrium: C (t) starts from zero at t = 
since an infinite temperature configuration is completely 
uncorrclated from one at equilibrium at temperature T . 
Afterwards, it increases and if the system reaches equilib- 
rium then it becomes equal to the typical correlation be- 
tween two equilibrium configurations, i.e. the self-overlap 
q computed in the previous section, or equivalently the 
long t — t' limit of C(t,t'), see top panel of Fig. 3. An 
asymptotic limit of C(t) below q instead indicates that 
the system is out of equilibrium even at long times. Actu- 
ally, this corresponds to the so called BIC test [16] , which 
has been introduced in recent numerical studies on sys- 
tems with pinned particles, where reaching equilibrium is 
a delicate issue. 

Guided by the knowledge of the dynamical behavior at 
c = 0, we expect a system quenched just slightly below the 
dynamical transition line to never equilibrate and to show 
the aging dynamics discovered in [T5]- Indeed, our nu- 
merical solution^ of (0| confirms this intuition (see Fig0]): 
C(t,t') strongly depends on the waiting time t' even for 
large time-scales, but a plot of the same curves in terms of 
the relevant adimensional variable 9 = t/t' shows a good 
collapse. At the same time, on time-scales of the order 
of one time translation invariance (TTI) and the fluctua- 
tion dissipation theorem (FDT) hold. On this basis and 



2 Note that the spherical parameter /i(t) is not related to the 
energy per spin in the same simple way found for the unconstrained 
p-spin spherical model. 

3 We used the efficient numerical integration method proposed 
by Kim and Latz in 1171 and developed in I18| . The main idea of 
this method consists in using an adaptive integration step to better 
describe the short and long time shape of the correlation function. 
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Fig. 3: Top panel: a successful BIC test showing that the 
asymptotic limit of the correlation function (as a function of 
t = t — t') and C(t) have the same large time value (T and c 
correspond to the point A3 of Fig[2]). Bottom panel: C(t) for 
several different concentrations of pinned spins at T = 0.65. 
The relaxation time rapidly grows approaching Cd f from the 
right. 



following the analysis |15j we simplify eqs. Q by sepa- 
rately investigating the TTI regime for which r = t — t' 
such that (t — t')/t — > and the aging regime for which 
(t — t')/t ~ 0(1) and where we use = t'/t. In the first 
regime the equation on C(t,t') is related to the one of 
R(t,t') by FDT and reads: 

d T C(r)+ (V(l) - Moo + (1 - C{r)) = (6) 



- -C(r) - p / dud u C(u)V'(C(T - «)) . 
P Jo 

Here fiaa is the large time value of the spherical parameter 
in the non-equilibrium regime. Taking the large r limit 
of ©, /ioo is found directly related to C(t — > oo) = n: 



PfJ-oo = 



1 



1-n 



/3 2 [V'(l)-V'(n)] 



(7) 



To obtain the value of one has to solve the aging regime 
since ^ depends on it, see (JSJ) - In order to do that, we 
first introduce the notation C(t, t') = r\C{ff) and R(t, t') = 
H{9). Since lim^oo C(t) = n, we have C(l) = 1. The 
t t' limit of the correlation function is denoted ro , hence 
C(0) = ro/ri. The two equations on C{6) and 72.(0), that 
we do not reproduce here, coincide under the hypothesis 
that C and 1Z obey to a generalized FDT relation 1Z(9) = 
xj3riC'(9), where x quantifies the FDT violation (x = 1 
means no violation). 

In summary, for a quench to a given point in the plane c— T 
five parameters should be determined to characterize the 
aging behavior: x, /ioo, Coo, ?T, and r$. The five equations 
needed for the solution are (0c]), (J5j) , (JT)) , and the two ones 
given by the equations on C(6) (or TV) evaluated in 9 = I 
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Fig. 4: Aging behavior of C(t,t') for a quench at the point A2 
of FigO The dotted line denotes the plateau value r%. 

and 6 = 0. The final set of equations includes ([7]) and 

/?Moo = 1+ (8a) 

+ p 2 [v'(i) + Coo V (Coo) - (i - x)nv'(n) - xr v'M] 

C oo = p 2 V(C oo )(l-r 1 (l-x)-xr ) (8b) 
l = /3 2 (l-n) 2 V"(ri) (8c) 
/Wo = P 2 [V'(r ) + r V'(l) + C 00 V'(C 7 oo)+ (8d) 
-(1 - x)nV'(ro) - (1 - x)r V'(n) - 2xr V(r )} . 

We have studied the solution of these equations varying 
the values of T and c. In order to avoid spurious results, 
we follow the solution starting from the dynamical transi- 
tion line, where one can show that n = qi, ro = qo = Coo, 
x = 1, and /loo = Mcq = 1 / P + /3V (1) ■ Starting from ca{T) 
and increasing c or decreasing T, we find that and Coo 
rise gently, x declines from one and r\ remains similar to 
qi (in general ro < Coo < ri). For a large enough value 
of c (or a small enough value of T) the solution ceases 
to exist. This leads to a second dynamical transition line 
Td (c), or equivalcntly a,{T), that meets the MCT dy- 
namical transition line at (ch,Th), see Fig. 2. Within the 
region delimited by the two dynamical transition lines the 
system never equilibrates and ages forever. In this regime 
the physical behavior is the one already understood for 
the unconstrained p-spin model. Actually, the free energy 
landscape of the randomly pinned p-spin spherical model 
is very similar to the one of the unconstrained system, ex- 
cept that there is one special state that is singled out and 
favored by the random pinning field pQ . After a quench be- 
tween Cd{T) and Cd f (T) the system ages because in its de- 
scent in the free-energy landscape, it is never able to drift 
lower than the so called marginal states. These "trap" the 
system: the closer it approaches them, the fewer are the 
unstable modes allowing it to escape. The existence of a 
preferred state correlated with the random pinning field 
is irrelevant since the system is not able to visit it in this 
regime. In order to be sure that this interpretation is cor- 
rect we also independently obtained ©, (|5aj) . (|8b|) . and 
(|8dj) using the replica formalism. In this context Coo is the 



overlap between marginal states and equilibrium states, 
ri and ro represents the self and the mutual overlap of 
marginal metastable states and x is the block parameter 
of the one step replica symmetry breaking matrix. 
Let us now discuss the nature of the transition at Cd s (T), 
which is different from the usual MCT one. In this case 
the existence of a preferred state makes a difference. In 
fact, beyond Cd f (T) the pinning field is so strong that the 
marginal states disappear and the system ends up in the 
equilibrium state favored by the field. However, contrary 
to what happens at c<j(T) the marginal states are not typ- 
ical equilibrium states, actually they are quite different 
from those. Thus, approaching Cd f {T) from the right the 
system ages for a certain time, since it is lost around what 
remains of the marginal states, and then at certain point 
it finds some escaping direction and relaxes to equilib- 
rium. Thus, as it happens for quenches close to the MCT 
transition line, aging eventually is interrupted but in a 
way which is very different. This is clearly shown in the 
lower panel of Figj3) C(t) first rises toward an interme- 
diate value, then it suddenly jumps to the equilibrium 
self-overlap. Concomitantly, C(t, t') stops to age and TTI 
and FDT start to hold. Approaching Cd f (T) the time 
to reach equilibration becomes longer and longer, possi- 
bly diverging as a power law in c — Cd f {T). In order to 
understand better the nature of this transition it is use- 
ful to introduce two different timescales: the equilibrium 
decorrclation time, r eq , which is the time it takes to the 
equilibrium correlation function C(r) to decrease to its 
asymptotic limit and the relaxation time, r re i, which is 
the time it takes to C{t) to reach the self-overlap equilib- 
rium value. Approaching Cd(T) from the left, r eq and T re i 
both diverge. Instead, approaching c<j, (T) from the right 
only r re i diverges: the system takes a long time to find 
the way to escape from (almost) marginal states but then, 
when this is found, the system relaxes to the equilibrium 
state which is very different from the marginal states and 
within which relaxation is fast. This kind of transition has 
been also found in some quantum spin disordered models 
|19| and dubbed first order dynamical phase transition. 
We think that physically it resembles more to a dynamic 
spinodal because it is related to the sudden disappearance 
of the marginal states. 

Discussion and Conclusion. The two main re- 
sults of the previous sections, based on mean field the- 
ory, concern the MCT phase diagram for randomly pinned 
glassy systems and the aging behavior for quenches in the 
c — T plane. Let us now discuss them, especially their ex- 
tension to finite dimensions and realistic systems. First, 
we would like to stress that although our results are ob- 
tained for an admittedly very abstract model, they are di- 
rect consequences of the evolution of the free-energy land- 
scape in the c — T plane. Thus, they are expected to be 
valid for all systems sharing a similar landscape and in 
particular, within RFOT theory, also for particle models 
in finite dimensions. As anticipated in the introduction 
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our MCT phase diagram differs from the one obtained 
for hard spheres by Krakoviack [TU] using a generalization 
of the MCT projection operators method. He finds that 
the transition line does not end in a critical point. In- 
stead, it goes on but the transition becomes continuous. 
On the contrary, after the end-point, we just find a Widom 
cross-over line where the decorrelation time is maximum 
but not infinitc[3- At this stage it is not clear what is 
the origin of this difference, besides the fact that our ap- 
proach only focuses on collective dynamical slowing down 
whereas Karkoviack's one somehow put together single- 
particle and many particles phenomena. However, this 
disagreement between the two approaches partially fades 
out when the results are transposed in finite dimensions 
since the MCT transition line becomes a cross-over any- 
way. It would be interesting to perform either the repli- 
cated HNC computation of [3U] or the cage expansion of 
[22[55] for randomly pinned systems to better clarify this 
issue and obtain quantitative results complementary to 
the ones in [TO] . 

Let us now focus on the dynamical behavior below Td(c). 
Since the MCT transition is avoided in finite dimensions, 
r re i only diverges at Tr-(c) and in a a super- Arrhenius 
way. Increasing c the divergence of the relaxation time 
becomes milder. Indeed at (c^T^) the exponent of the 
generalized Vogel-Fulcher-Tamman law becomes equal to 
1/6 within mean-field theory (and 1 within the Migdal- 
Kadanoff RG procedure of pQ) while it is equal to 2 at 
c = pp. Thus, the system becomes stronger and, as dis- 
cussed previously, less dynamically heterogeneous. This 
is in agreement with the results of numerical simulations 
[5]. We now consider the aging dynamics. For an un- 
constrained system, quenches above Tr-(O) lead to inter- 
rupted aging whereas quenches below it lead to aging that 
goes on forever. Instead, for constrained ones, we expect 
interrupted aging also below Tk{c) (in particular above 
Td f (c)). The reason relies in the form of the free energy 
landscape: after the quench the system is trapped and ages 
on high free-energy states in way very similar to what hap- 
pens above Tr-(c). However, below Tk(c), the state corre- 
lated with the pinning field is thermodynamically favored 
with respect to all the others even taking into account 
their multiplicity. Thus, when the waiting time becomes 
of the order of the nucleation time of the favored state the 
system stops aging and suddenljjf] reaches equilibriurn^l, 
accordingly the decorrelation time within the equilibrium 
state is very different from the relaxation time, the former 
is expected to remain small whereas the latter diverges 
approaching Tff(c) from below. This is relevant for the so 

4 Single particle Lorentz gas physics and percolation phenomena 
may play a role at high temperature, where higher concentrations 
of pinned particles are required to reach the interesting region of 
the phase diagram. These could make the observation of collective 
dynamical slowing down more difficult. 

5 Actually the nucleation process can be slow because it takes 
place in a random environment. 

"Remarkably, a similar phenomenon plays an important role in 
optimization problems such as compressed sensing 1231 . 



called BIC test: it means that in this regime starting from 
the pinned configurations the system is automatically at 
equilibrium and decorrelates rapidly 0. Thus one can eas- 
ily obtain equilibrated averages even though the time it 
takes to get a successful BIC test could be very largc^. 
Numerical simulations and even experiments in colloids 
using optical trapping techniques can be performed to 
test our predictions. In particular, it would be very in- 
teresting to observe the nucleation-interrupted aging be- 
low Tr-(c) and the characteristic dynamical features taking 
place near the spinodal point (cfc,T^). 

* * * 
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